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ABSTRACT 

Intrinsic alignments of galaxies can mimic to an extent the effects of shear caused by weak 
gravitational lensing. Previous studies have shown that for shallow surveys with median 
redshifts z m ^0.1, the intrinsic alignment dominates the lensing signal. For deep surveys 
with z m ~ 1, intrinsic alignments are believed to be a significant contaminant of the lensing 
signal, preventing high-precision measurements of the matter power spectrum. In this paper 
we show how distance information, either spectroscopic or photometric redshifts, can be 
used to down-weight nearby pairs in an optimised way, to reduce the errors in the shear 
signal arising from intrinsic alignments. Provided a conservatively large intrinsic alignment 
is assumed, the optimised weights will essentially remove all traces of contamination. For 
the Sloan spectroscopic galaxy sample, residual shot noise continues to render it unsuitable 
for weak lensing studies. However, a dramatic improvement for the slightly deeper Sloan 
photometric survey is found, whereby the intrinsic contribution, at angular scales greater than 
1 arcminute, is reduced from about 80 times the lensing signal to a 10% effect. For deeper 
surveys such as the COMBO-17 survey with z m ~ 0.6, the optimisation reduces the error 
from a largely systematic ~ 220% error at small angular scales to a much smaller and largely 
statistical error of only 17% of the expected lensing signal. We therefore propose that future 
weak lensing surveys be accompanied by the acquisition of photometric redshifts, in order to 
remove fully the unknown intrinsic alignment errors from weak lensing detections. 

Key words: cosmology: observations - gravitational lensing - large scale structure, galaxies: 
formation 



1 INTRODUCTION 

Weak gravitational lensing by intervening large-scale structure in- 
troduces a coherent distortion in faint galaxy images. Several inde- 
pendent surveys have measured this 'cosmic shear' effect and are 
now able estimate the bias parameter, b, (Hoekstra, Yee & Gladders 
2001; Hoekstra et al. 2002a) and set joint constraints on the values 
of the matter density parameter Q. m and the amplitude of the matter 
power spectrum, a% (Bacon et al. 2002; Brown et al. 2002a; Hoek- 
stra, Yee & Gladders 2002; R.Maoli et al. 2001; Rhodes, Refregier 
& Groth 2001; Van Waerbeke et al. 2001). 

A potential limitation of this technique is the intrinsic correla- 
tion of the ellipticities of nearby galaxies which can result from grav- 
itational interactions during galaxy formation. This 'tidal torquing' 
leads to a net alignment of nearby galaxies which could produce a 
spurious signal similar to that induced by weak gravitational lens- 
ing. A number of studies have investigated the amplitude of this 
intrinsic alignment effect with some broad agreement, (Heavens, 
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Refregier & Heymans 2000, hereafter HRH, Brown et al. 2002b; 
Catelan, Kamionkowski & Blandford 2001; Crittenden et al. 2001; 
Croft & Metzler 2000; Hui & Zhang 2002; Lee & Pen 2001; Jing 
2002; Porciani, Dekel & Hoffman 2002), but it is fair to say that an 
accurate estimate of this effect eludes us. For low-redshift surveys, 
for example SuperCOSMOS and the Sloan spectroscopic galaxy 
sample, where the median redshift is z m ~ 0.1, the correlation of 
ellipticities due to intrinsic alignment is far greater than the expected 
lensing signal. As the survey depth increases, galaxies at fixed angu- 
lar separation are distributed over wider ranges of redshifts, leaving 
only a small proportion of the galaxy pairs close enough for tidal 
interactions to correlate the ellipticities. This reduction in the to- 
tal angular intrinsic alignment signal, combined with an increasing 
line-of-sight distance boosting the lensing signal, leaves a smaller 
intrinsic alignment contamination for deeper surveys. The majority 
of studies find that for surveys with depths z m ~ 1 the estimated 
contamination contributes to up to ~ 10% of the lensing signal, al- 
though Jing (2002) has argued that the contamination could be much 
higher. As cosmic shear measurements become more accurate these 
low levels of contamination cannot be considered negligible. 

If an accurate estimate of the intrinsic ellipticity correlation 
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strength did exist, then it would be possible to subtract it, leav- 
ing the correlation induced purely by lensing. In the absence of a 
good estimate for the intrinsic correlation function, it is obvious that 
one can improve upon a straightforward ellipticity correlation by 
down-weighting galaxy pairs close in redshift and angular sky sep- 
aration. This can be done at the expense of increasing the shot noise 
contribution from the distribution of individual galaxy ellipticities. 

In this work, we deduce the optimal pair weighting for an arbi- 
trary survey, given accurate information about the galaxy distances. 
We use these results on the Sloan spectroscopic survey design to 
show that the optimal weighting scheme can largely remove the con- 
tamination by intrinsic alignments leaving almost pure shot noise. 
The shallow depth of this spectroscopic sample still prevents its use 
for weak lensing studies as the remaining shot noise exceeds the 
expected weak lensing signal. 

For multi-colour surveys with photometric redshift informa- 
tion we propose a semi-optimised procedure of excluding pairs 
which are likely to be close in three dimensions. Comparison with 
the optimal method shows this can be very effective. We apply this 
technique to the Sloan photometric survey (SDSS) with z m ~ 0.2, 
and three deeper surveys, the Red-Sequence Cluster Survey (RCS) 
with z m ~ 0.56, (Hoekstra et al. 2002b), a sample of the COMBO- 
17 survey with z m ~ 0.6, (Wolf, Meisenheimer & Roser 2001) 
and the Oxford Dartmouth Thirty survey (ODT) with z m ~ 1.0, 
(Allen, Dalton & Moustakas 2002). We find that even with fairly 
inaccurate photometric redshifts it is possible to reduce the contam- 
ination from intrinsic alignments significantly. For the photometric 
SDSS, the improvement is dramatic, enabling it to be used as a 
weak lensing survey. For the deeper surveys we show that using 
this weighting scheme, intrinsic alignment contamination can be 
reduced by several orders of magnitude. 

The effect of the non-uniform weighting on the weak lensing 
shear correlation function is calculated. We find that galaxy pair 
weighting slightly reduces the lensing signal, the reduction depen- 
dent on the photometric redshift accuracy and the survey depth. 

This paper is organised as follows. In §^we briefly describe re- 
lated weak lensing theory. In §^ we derive the optimal pair weighting 
for a spectroscopic survey. We discuss the results of pair weighting 
for the Sloan spectroscopic sample in comparing HRH and Jing 
intrinsic alignment models. In §n we propose an alternative scheme 
for multi-colour surveys with photometric redshift information and 
present the results in These results are compared with the ex- 
pected weighted lensing signal derived in the Appendix. In §^ we 
discuss the results and the implications for the design of future weak 
lensing surveys. 



orientated galaxies in the absence of weak lensing, a correlation of 
galaxy ellipticities would yield a zero result at all angular scales. The 
presence of weak lensing introduces coherent shearing distortions, 
producing ellipticity or shear correlations which are related to the 
nonlinear mass power spectrum Pg via 
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where Jo is the zeroth Bessel function of the first order and P K (k) 
is the convergence power spectrum at wave number k, 
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}k (w) is the comoving angular diameter distance out to a comoving 
radial geodesic distance w, and is defined in the Appendix; a(w) is 
the dimensionless scale factor, Ho is the Hubble parameter and Q m 
the matter density parameter. The second argument of Ps allows for 
time-evolution of the power spectrum. g(w) is a weighting function 
locating the lensed sources, 



g{w)= / dw 4>{w) J ; ' . 



(4) 



<f>(w(z)) is the redshift distribution or selection function, and wh 
is the horizon distance, (Bartelmann & Schneider 2001). 

In addition to shear correlations produced by lensing, several 
groups have found evidence for an intrinsic alignment component. 
This means that to some degree, galaxies are not randomly orien- 
tated. In the absence of intrinsic correlations, the shear correlation 
function (77* ) can be estimated from a catalogue of galaxy shapes 
through estimates of the ellipticity correlation function 
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where e a is the ellipticity of galaxy a. In practice the pair sum is 
taken over galaxies in a small range A9. 

Typically the weight assigned to each pair of galaxies, W a b, is 
the product of the uncertainty in each galaxy shape measurement. In 
the next section we derive an additional optimal pair weight which, 
in effect, down-weights galaxy pairs close enough to contribute to 
a contaminating intrinsic alignment signal, enabling the elliptic- 
ity correlation estimator to be directly related to the lensing shear 
correlation function. 



2 THEORY 

In the weak lensing limit the ellipticity of a galaxy, (ei, ei), is an 
unbiased estimate of the shear 7, the stretching and compression 
of the galaxy image caused by the gravitational lensing of its light. 
Ellipticity is defined by approximating each galaxy as an ellipse 
with axial ratio /3, at position angle <j>, measured counterclockwise 
from the x-axis. 
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The complex shear 7 is related to the average galaxy ellipticity 
by 7 ~ (e)/2 = (ei + ie%)/2. A useful statistic which can be 
used to constrain the matter power spectrum is the shear correlation 
function, (7*y*)e, for galaxies separated by angle 9. Forrandomly- 



3 OPTIMAL PAIR WEIGHTING FOR SPECTROSCOPIC 
SURVEYS 

The ellipticity correlation function, for pairs separated by an angle 
9 is given by equation ^[ This estimate includes an uncertain con- 
tribution from intrinsic alignments, I a b. The estimate of the lensing 
shear correlation function £[77*; 9] is therefore 
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The amplitude of the intrinsic correlation is controlled by the co- 
moving distance between the galaxy pair. It is only significant for 
galaxies closer than a few tens of Mpc at most, in which case the 
comoving separation is given by 
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In principle, the intrinsic ellipticity correlation depends on the ori- 
entation of the galaxy pair in three dimensions, and also on the 
redshifts of the galaxies. However, studies of the intrinsic effect 
average over the orientation, and there seems to be little evolution 
with redshift. Thus we assume we have an estimate of an average 
correlation function which depends only on the separation 7? of the 
two galaxies: 



r)i(R) = (e(x)e*(x + R))i 



(8) 



where the average is over galaxy positions x, assumed to be known 
precisely. In practice, peculiar velocity distortions and spectroscopic 
redshift uncertainties lead to errors of order Az ~ 0.001, which we 
will neglect. This is justified because, as we show in section 6, sur- 
veys with photometric redshift errors much larger than this can have 
intrinsic signals removed virtually as effectively as the idealised ex- 
ample treated here, rji (R) can be calculated approximately, from 
numerical simulations (HRH, Croft & Metzler 2000; Jing 2002) 
or analytic approximation (Catelan, Kamionkowski & Blandford 
2001; Crittenden et al. 2001; Hui & Zhang 2002; Lee & Pen 2001; 
Porciani, Dekel & Hoffman 2002). The published results do differ 
by more than an order of magnitude, and this uncertainty has to be 
considered when choosing which rji (R) model to use in the follow- 
ing pair weighting scheme. The preferred strategy when applying 
this technique is to assume the largest theoretical model for rji (R) 
to ensure that all feasible contamination from intrinsic alignments 
is removed. This is discussed further in section 6. 

To calculate the angular intrinsic contribution to the weak lens- 
ing signal, rji (R) is averaged over galaxy pairs separated by 9 with 
a redshift distribution characteristic of the depth of the survey. This 
can be compared to the signal measured from low-redshift surveys 
such as the SuperCOSMOS survey (Brown et al. 2002b), where el- 
lipticity correlations are dominated by intrinsic alignments. These 
studies show very rough agreement of the amplitude of the effect, 
but the problem is sufficiently difficult that a definitive accurate 
study has yet to emerge. To quantify our uncertainty on the intrin- 
sic correlations, we assume the intrinsic alignment prediction has a 
fractional error / and thus a variance 



°1A = f{Iab) 2 
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where {I a b) is a weighted average, and we will conservatively take 
/ ~ 1 in what follows. The error in the estimator of the lensing 
correlation function, equation ^, has two main sources: uncertainty 
in the intrinsic correction, and shot noise from the intrinsic ellipticity 
distribution. If the variance of Real(e a ej*) is fTp air , then the shot noise 



(E ab w ab ) 

In the limit of an infinite number of galaxies, <tsn tends to zero. 

Considering the shot noise error and intrinsic alignment error 
only, the total error on the lensing correlation function is then given 
by 
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We minimise this variance subject to ^2 ab W ab = constant. The 
weights for a galaxy pair p — {a, b} then satisfy 
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where A' is a Lagrange multiplier. If we define a matrix 



where 
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then W p satisfies 
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where A is an unimportant constant which we will set later by 
choosing the weights of uncorrelated pairs ( J p =0) to be unity. U p 
is a 'unit' vector consisting entirely of ones. The solution is 



W= (I + M) _1 U. 



(16) 



Because of the form of M, the inverse can be computed, by 
expanding (I + M) _1 = I - M + M 2 - . . ., with elements 
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The weight is then obtained from equation (jlfj), equivalent to sum- 
ming (Il8|) over pairs q: 



W„ = 1 - 



1 +A2 



(19) 



where Ai = ^ J r and we have chosen A = 1. This is the optimal 
weighting scheme if distances to the galaxies are accurately known. 



3.1 Reduction in variance 

The variance in the lensing signal is given by equation (pi]), which 
may be written 
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For comparison, with equal weighting of galaxy pairs, the vari- 
ance is 
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For optimal weighting, 
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The ratio of the optimal variance to the equal weighting variance is 
shown in Fig. 1 as a function of Ai = Ai / v~N and A2, and is given 
by 



Ratio 



(I+A2) 



(1 + A 2 )(1+A 2 -A 2 )- 



(23) 



Note that Aj < A2, which can be shown by considering the positive 
quantity J2 p ( J p - ( J p)Y- 

The sums over galaxy pairs may be estimated by integrals 
over the mean number density of the survey, assuming that we 
have accurate comoving galaxy distances which could come from 
a spectroscopic redshift survey such as the Sloan spectroscopic 
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Variance improvement 
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Figure 1. Ratio of optimal variance to the equal weighting variance. 

galaxy sample. Redshifts can be converted into distance coordinates 
w using the approximate formulae in Pen (1999). 
In order to set the weights, we need to compute 
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where the selection function (j) w (w) is the number of objects at 
distance w included in the survey, and £ is the galaxy two-point 
correlation function. <j> w is related to the redshift selection function 
(j>z by <j> w dw — <j) z dz and we use the redshift distribution of the 
form. 



4> z {z) oc z exp 



(25) 



where the median redshift of the survey is z m ~ IAzq. This is 
normalised by J dzcj> z (z) — n s i r , where n str is the number of 
galaxies per steradian in the survey. Note that for n = 0we calculate 
N, the expected number of pairs at separation between 8 and 8+A8. 
For (Tp a i r , we assume that the dispersion in ellipticity is a e ~ 0.2, 
(Hudson et al. 1998), giving crp ail - ~ 0.04. 

Typically, for angular scales less than 8 ~ 15 arcminutes, 
we find that Aj <C A2. The ratio, equation [23[ is then simply 
Aj~ 2 , and the optimally-weighted error is almost pure shot-noise 
a\ ~ (Tpair/A^. Changing A8 predominantly alters and hence the 
amplitude of the shot noise component. 

Note that the weights are not constrained to be positive and 
in extreme cases, if J p is large, can be negative. The reason for 
this is that the contribution to the error from intrinsic alignments is 
assumed to be proportional to the intrinsic alignment signal itself. 
If this is the dominant error term, it is possible that one can have a 
net gain by reducing this term through giving some pairs negative 
weight. This behaviour is, however, rare. 



4 RESULTS FOR THE SLOAN SPECTROSCOPIC 
SAMPLE: COMPARISON OF HRH AND JING 
INTRINSIC ALIGNMENT MODELS 

Covering n steradian the Sloan Digital Sky Survey, SDSS, is a wide, 
shallow survey. The spectroscopic galaxy sample has a median 
redshift of 0. 1 with redshift measurements known to an accuracy 
of Az = 0.0001 (Stoughton & SSDS Team 2002). Figure § shows 
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Figure 2. Reduction in the error from intrinsic alignments and shot noise 
for the Sloan spectroscopic sample. The upper curves show the error in 
the shear correlation function for two intrinsic alignment models, assuming 
all galaxy pairs are weighted equally. Note that these curves essentially 
combine a statistical error (shot noise) with a systematic contamination 
(intrinsic alignments). Triangles show the Jing model and circles HRH's 
spiral model. The error with optimal weighting is shown as rilled circles, 
and has reduced both intrinsic alignment model signals close to shot noise 
(shown dashed). Note that the increase in shot noise is negligible due to the 
small proportion of downweighted pairs in the optimal weighting scheme. 
The expected weak gravitational lensing correlation function for a ACDM 
model, (dotted), is still dominated by shot noise for this shallow z m =0.1 
sample. 



the reduced intrinsic alignment error calculated using the optimal 
weighting scheme, compared to the error introduced to the lensing 
signal by an unweighted HRH intrinsic alignment signal and an 
unweighted Jing intrinsic alignment signal at a redshift of 0. 1 . The 
minimised signal for both intrinsic alignment models is close to 
pure shot noise with residuals of less than 1%. 

The published differences between the HRH and Jing models 
come from from the higher resolution N-body simulation used by 
Jing, and differences in the choice of £ ( R) . Jing finds that the spatial 
ellipticity correlation of the halos, r\i (R), derived for the same min- 
imum dark matter halo mass, is higher by a factor of 2, compared to 
the elliptical HRH model. This model assumes that the ellipticity of 
the galaxy is the same as that of the parent halo. It is sensitive to the 
accuracy of the halo shape measurement and is therefore better de- 
termined in the higher resolution simulations. A second difference, 
in the angular intrinsic ellipticity correlation, arises since Jing uses a 
larger two-point correlation function, measured from the simulation 
itself, whereas HRH assumed £(R) = (R / Sft-'Mpc) -1 ' 8 . HRH 
also used a model for spiral galaxies, where the galaxy is modelled 
as a thin disk, perpendicular to the angular momentum vector of the 
parent halo. The spiral model has significantly lower correlations at 
spatial separations R a i, > 5h~ Mpc. 

Except where otherwise stated, we use the HRH spiral model 
throughout this paper. r\i (R), summed over e\ and e%, is fitted by an 
exponential: r]i(R) = 0.012exp(-i?/[1.5/i -1 Mpc]). We should 
point out that fitting an exponential model could underestimate 
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Survey 


size str 




A z 


Sloan spectroscopic sample 


7T 


0.1 


0.0001 


SDSS photometric 


IT 


0.2 


0.025 


COMBO-17 


0.0005 


0.6 


0.03 


RCS BVRz' 


0.015 


0.56 


0.3 


ODT UBVRIz' 


0.005 


1.0 


0.2 



Table 1. Survey parameters for SDSS, COMBO-17, RCS and ODT 

rji(R) on large scales R > 10h~'Mpc, where the signal-to-noise 
is low. This would then lead to an underestimation of the intrinsic 
signal on angular scales 6 > 10'. If this were the case, however, the 
results would be inconsistent with the lack of detection of a B mode 
signal on these scales found by Hoekstra, Yee & Gladders (2002), 
discussed further in section 6. Figure |] shows the contamination 
expected from this exponential HRH model and the larger signal 
found by Jing (2002), assuming the same minimum halo mass as 
HRH (6.9 x 1O U M0). We have, however, used the power-law two- 
point correlation function to convert the Jing rij(R) fitting formula 
into an angular signal. Increasing the minimum halo mass increases 
the unweighted Jing signal, but does not greatly affect the minimised 
signal leaving residuals at less than 3% of the shot noise. 

As the Sloan spectroscopic sample is so shallow, the expected 
lensing signal is negligible preventing this spectroscopic sample 
from being used as a weak lensing data set. It is however an inter- 
esting test of the method showing that the optimal weighting scheme 
will almost completely remove both the HRH and the stronger Jing 
model for intrinsic alignments. We will therefore proceed with the 
HRH spiral model for intrinsic alignment, as the method has been 
shown to remove successfully any amplitude of intrinsic alignment 
provided the signal is significant only for galaxy pairs closer than a 
few tens of Mpc. 



5 USING PHOTOMETRIC REDSHIFTS 

The method detailed in section [5] is only feasible for spectroscopic 
surveys with essentially known galaxy distances. More typically we 
have good estimates of galaxy distances from photometric redshifts 
with associated errors. These errors are much larger than the scale 
over which the intrinsic alignments are correlated. We must there- 
fore define a new weighting scheme for galaxy pairs dependent on 
their estimated redshifts alone. Those weights are then transformed 
into an average weight assigned to a pair at true redshift separation 
to calculate the residual angular intrinsic contribution. 

Using the minimum lensing error derived in section |as a 
benchmark we propose a simpler approach to derive the weight- 
ings for galaxies at estimated distances. Applying the accuracy of 
the Sloan spectroscopic sample, A z = 0.0001, to this alternative 
method we come very close to reproducing the previously derived 
minimum. This shows that our alternative method, in the limit of 
accurate photometric redshift information, is close to the optimal 
result. As the errors in the distance measurements increase, our 
ability to down-weight close galaxy pairs correctly decreases, and 
the residual intrinsic contamination increases. 

For a pair of galaxies with estimated redshift z a and z b we 
assign a zero weight if \z a — z b \ < aA z and a weight of one 
otherwise. We choose a to minimise the total error on the shear 
correlation function. The optimum value depends on angular sep- 
aration, and survey details. Clearly, this is not the most general 
weighting scheme, but we will see later that, with good photometric 
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Figure 3. Residual error expected after removing pairs within aA z of each 
other for the ODT survey. A z is the root mean square photometric redshift 
eiTor. 



redshifts, this simple procedure does almost as well as the theo- 
retical optimum where all galaxy distances are known. Assuming 
that the estimated redshifts are distributed normally about the true 
redshifts with Gaussian widths given by A z , we find that the average 
weight for a galaxy pair at true redshifts z a and Zb is given by 



{Wab) = 1 ^ 

where 

Za Z a 



1 



20F 



dxe. x [erf(j/ -+- a:) — erf(u + x)] (26) 



x = 



\/2A z 



- Zb + aA z 
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z a — Zb — aA z 
V2A Z 



(27) 



The weighted intrinsic component can then be evaluated using the 
integral version of equation [ill setting the value of a by minimising 
the total error contribution from the intrinsic alignments, ctia and 
the shot noise, <tsn- 

_ // dz a dz b (f) z (Za)<t> z (zb) [l + £{R a b)] {W a b)Vl (Rab) 
J dz a dz b (j> z (Za)<t) z (z b ) [l + £{R a b)] {W a b} 

> T <\ = ' T n;;i, j / dz a dz b (j) z (Za)(t> z {zb) [l + £ {Rab)] {W ab )^j (29) 

As a increases the probability of removing all close galaxy pairs in- 
creases and the intrinsic contribution decreases to zero. This means 
that the total galaxy count decreases and the shot noise increases, 
hence there is an optimum value of a. As 8 increases, shot noise 
begins to dominate the intrinsic alignment error, and a progressively 
lower value of a is favoured. This is shown in figure ^| where the 
reduced intrinsic contamination is plotted for varying values of a 
for a low accuracy photometric redshift survey, for example the 
ODT. 

The optimum values of a are dependent on survey size, depth, 
photometric redshift accuracy and the choice of rji(R). 
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6 APPLICATION OF SEMI-OPTIMISED WEIGHTING 
SCHEME TO MULTI-COLOUR SURVEY DESIGNS 

We consider four multi-colour surveys with varying depth, angu- 
lar size and photometric redshift accuracy: Sloan photometric sur- 
vey, SDSS, COMBO-17, the RCS survey and the ODT survey. 
Table |l| details the relevant survey parameters. The uncertainty 
in the redshift depends mainly on the number of colours known 
for each object. Photometric redshift errors for SDSS assumes a 
neural-network-based photometric redshift estimator (Tagliaferri 
et al. 2002; Firth, Lahav & Somerville 2002). For the RCS and 
ODT we assume the use of hyper-z (Bolzonella, Miralles & Pello 
2000). The full COMBO-17 survey has imaging to R ~ 25.5, but 
colour information in 17 bands (producing very accurate photomet- 
ric redshift estimates) down to a limiting magnitude of R ~ 24. 
It is this magnitude-limited sample of the full COMBO-17 weak 
lensing survey which we will consider. 

Figure ^ shows the correlation signal we expect to find from 
HRH-derived intrinsic alignments and the best reduction that can be 
achieved with this method, with photometric redshift information at 
the accuracy appointed for each survey. These have been calculated 
using the values of a(8) shown in figure ^. The dashed line in 
Figure ^ shows the optimal improvement in the error obtainable 
if accurate distances were known, using the method detailed in 
section ^. This acts as a useful benchmark to see how close the 
semi-optimal method for multi-colour surveys gets to the ideal. 
The intrinsic alignment signals ca n be compared to the expected 
weighted lensing signal, equation \13, for each survey. This is 



SDSS z„ = 0.2 



RCS survey z m = 0.56 




derived from a nonlinear CDM mass power spectrum with Q. m = 
0.3, Q. A = 0.7, T = 0.21, erg = 0.9, calculated using the 'halo 
moder fitting formula (Smith et al. 2002). For the pair weighting 
proposed, this signal is slightly lower than the unweighted lensing 
signal (equation □), (see Appendix). 

The true amplitude of the intrinsic alignment signal is rather 
uncertain; the correlation signal from intrinsic alignments could 
move up, as argued by Jing (2002), or down, as argued, for example, 
by Crittenden et al. (2001), and by up to an order of magnitude. 
Figure ^ also shows the optimum a(8) if the intrinsic alignment 
signal is as high as found by Jing (2002), for the ODT survey. Due 
to the presence of intrinsic correlations at larger galaxy pair spatial 
separations, it is necessary to remove more pairs at larger angular 
scales. At the expense of a higher residual shot noise, this more-or- 
less guarantees removal of any plausible intrinsic alignment signal, 
leaving errors at less than 3% of the weighted ODT lensing signal. 
This is probably the preferred strategy until the intrinsic signal is 
better determined, but it should be noted that an intrinsic alignment 
signal of the amplitude found by Jing (2002) is incompatible with 
current weak lensing detections, (Van Waerbeke, Mellier & Tereno 
2002). 

Using an HRH model for intrinsic alignments we find that the 
SDSS is dominated by intrinsic alignments and that the weak lensing 
signal from the two middle-depth surveys, RCS and COMBO-17, 
will be strongly contaminated at angular scales less than 10 arcmin- 
utes. Interestingly Hoekstra, Yee & Gladders (2002) show an E/B 
mode decomposition for the RCS survey which shows a significant 
B mode at angular scales less than 10 arcminutes. The presence of 
a B mode can be due, on small angular scales, to source cluster- 
ing, (Schneider, Van Waerbeke & Mellier 2002), or it can imply 
systematics within the data and/or the presence of intrinsic galaxy 
alignments (Crittenden et al. 2002). 

The application of the weighting scheme produces encouraging 
results, the most startling of which comes from the SDSS. Due to the 
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Figure 4. Reduction in the error from intrinsic alignments and shot noise for 
photometric SDSS, RCS, COMBO-17 and ODT. The semi-optimal weight- 
ing, (filled circles), has reduced the unweighted HRH intrinsic alignment 
error (circles) to well below the expected amplitude of the weighted weak 
lensing shear signal, (dotted). The effect of semi-optimal weighting can be 
compared to the optimally weighted error, (dashed), attainable with spec- 
troscopic redshifts. 



wide sky coverage and accurate photometric redshift information 
it is possible to reduce the intrinsic alignment signal from well 
above, to well below the expected lensing signal. This enables, in 
principle, the extraction of a lensing signal at angular scales, 9 > 1', 
with errors less than 10%. We see similar reductions for the other 
three surveys, with the success of the reduction dependent on the 
accuracy of the redshift estimates. 



7 CONCLUSIONS AND IMPLICATIONS FOR FUTURE 
WEAK LENSING SURVEYS 

In this paper we have shown how distance information can be used 
to reduce the contamination of the weak gravitational lensing shear 
signal by the intrinsic alignments of galaxies. Our principal find- 
ing is that the level of intrinsic ellipticity correlation in the shear 
correlation function can be reduced by up to several orders of magni- 
tude, by down-weighting appropriately the contribution from pairs 
of galaxies which are likely to be close in three dimensions. 

For spectroscopic surveys we have derived an optimal galaxy 
pair weighting which reduces the contamination to a negligible 
component, leaving almost pure shot noise. Application of the tech- 
nique to the relatively shallow Sloan spectroscopic galaxy sample 
reduces the error by up to two orders of magnitude, but still leaves 
the lensing signal undetectable, dominated now by shot noise. 

For multi-colour surveys with photometric redshift estimates 
we have proposed a partially-optimised method, which removes 
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Figure 5. Optimum values of a for the four different surveys, assuming 
the HRH spiral model for intrinsic alignments. The dashed line shows the 
optimal values of a for the ODT survey, assuming the higher intrinsic signal 
found by Jing (2002). At the expense of marginally-higher residual shot 
noise, this conservative pruning of the pairs virtually guarantees almost 
complete removal of the intrinsic alignment contamination. 

pairs with close photometric redshifts from the computation of the 
shear correlation function. We find that with accurate photomet- 
ric redshifts this simple method is almost as effective as the fully 
optimised method. We show that, even with relatively crude photo- 
metric redshift estimates, the contamination by intrinsic alignments 
can be significantly reduced. 

Similar conclusions were independently drawn by King & 
Schneider (2002) who simultaneously proposed a weighting scheme 
based on photometric redshifts. Our technique however requires the 
assumed knowledge of a model for the intrinsic alignments. Pro- 
vided one is conservative, assuming the largest feasible model, this 
weighting scheme will reduce the true intrinsic alignment contam- 
ination to lensing correlation signals to a negligible level. 

We have applied the method to four multi-colour survey de- 
signs. The shallow photometric SDSS survey with z m ~ 0.2, shows 
the most dramatic of improvements. The intrinsic alignment signal 
is expected to dominate completely the weak lensing shear signal 
from a survey of this depth, but with judicious removal of pairs 
within ~ 0. 14 in estimated redshift, the error from shot noise and 
intrinsic alignments is reduced to only ~ 10% of the expected 
shear signal. This opens up the possibility of using the SDSS for 
future cosmic shear studies. With deeper surveys, such as the RCS, 
COMBO- 17 and the ODT survey, intrinsic alignments in the weak 
lensing shear signal are significantly reduced. With the highly accu- 
rate photometric redshifts of COMBO-17, the reduction is close to 
optimal, reducing the largely systematic ~ 220% error at small an- 
gular scales to a much smaller and largely statistical error of ~ 17% 
of the expected lensing shear correlation. This limiting error is pre- 
dominantly shot noise which will decrease as the survey size grows. 
The ODT and RCS have less accurate distance estimates but even 
with current photometric redshift estimates, which undoubtedly will 
improve with time, the reduction is quite significant. For the wide 



RCS, the 130% contamination at angular scales 9 ~ 2', is reduced 
to an error ~ 20%. For deep multi-colour surveys such as the ODT, 
the intrinsic alignment signal is almost completely removed, leaving 
noise at a level of < 2% of the lensing signal. 

The aim for weak lensing surveys to date has been to go as 
deep and as wide as possible. With an increase in depth comes an 
increase in the expected lensing signal and for purely geometrical 
reasons we also see a decrease in intrinsic alignment correlations. 
To reduce shot noise, the error introduced by the intrinsic elliptic- 
ity distribution of galaxies, surveys can go wide and/or deep. This 
increases the number of lensed sources, in the limit of an infinite 
number of galaxies the shot noise goes to zero. One additional 
source of error so far unmentioned is cosmic variance which, for 
small area surveys, can dominate on scales larger than a few ar- 
cminutes, (Schneider et al. 2002). This can be reduced by sampling 
different areas of the sky and combining to produce a wide field 
survey. Note that downweighting close galaxy pairs affects this ran- 
dom error little, whilst effectively removing the systematic errors 
from intrinsic contamination. 

As the size of future lensing surveys increase, control of sys- 
tematic errors becomes increasingly important. Intrinsic alignment 
contamination is a potentially large source of error, even if it is 
at a level much lower than assumed here. This paper has shown 
that with the addition of accurate photometric redshift information, 
the presence of intrinsic alignments can be effectively removed. 
Hoekstra et al. (2002b) have shown that it is possible to extract a 
low weak lensing signal from a survey with a median redshift of 
0.56 and argue that using this shallower survey enables more ac- 
curate star-galaxy separation and provides a fairly well determined 
redshift distribution for the sources, another uncertainty in weak 
lensing measurements to date. The addition of depth information to 
weak lensing studies also opens new possibilities for the applica- 
tion of lensing tomography (Hu 2002; Hu 1999), and reconstruction 
of the full 3D cosmological mass distribution (Taylor 2001; Hu & 
Keeton 2002), but here we have concentrated only on its benefits 
in reducing intrinsic alignment contamination. In view of the dra- 
matic improvements it offers lensing signal detection we therefore 
propose that for future weak lensing surveys, emphasis should be 
placed on acquiring multi-colour data in a wide area. In the case of 
restricted telescope time, moderate depth surveys, z m ~ 0.6, may 
then become a more attractive alternative to ultra-deep observations, 
z m > 1.0, that can often be limited by seeing. 
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APPENDIX A: THEORETICAL LENSING SIGNAL FOR WEIGHTED ELLIPTICITIES 

The effective convergence is (Bartelmann & Schneider (200f ), equation 6.18) 

= A [°° d Wg { W )f K { W ) 5[fK ^ M (Al) 
Jo a ( w > 

where A = 3H$Q. m / (2c 2 ), g( w) is the weighting function given in equation ^, and fx(w) is the geometry-dependent comoving angular 
diameter distance, dependent on the curvature K , 

( K-^ 2 sm(K'^w) (K>0) 
f K (w) = I w (K = 0) . (A2) 

[ (-/sT) _1/3 sinht(-- K ') I/2 * ] ( K <°) 

Note that w, the comoving radial geodesic distance, plays two roles, both as a third spatial coordinate, and as a time evolution label. 

The effective 2D convergence is an integral of the 3D convergence, and it is this which we will weight, since each galaxy gives an 
estimate of the shear, or convergence, in the weak lensing limit. 

/>oo 

«eff( 9 ) = / dw 4>(w)Kes( 9 , w) . (A3) 
Jo 

We will weight 3D convergence estimates with a weight function H(w, w' , | 6 ), to give the weighted correlation function 

{Kcb( )«eff( 9 )) = — I dw dw' H(w,w' ,| 9 - 9 \)(Kes( 9 ,w)n* ff ( 9 ,to'))<f>(w)(j>(w') 
y Jo Jo 

= b r dw r dw' 4>(wu(w')Hiw, W ',\t -Yd r d Wl fK{m \ fK } w \ m) 

Jo Jo Jo f K {w)a(wi) 

, fK(w 2 )f K (w' -W 2 ) . , , , , 

dw 2 . ,. . — : (o(fic[wi) 9 , Wild (Jk(w2) 9 ,w 2 )) (A4) 

j K (w')a(w 2 ) 

where Q = dwdw' H(w,w', 9)<j>{w)(j){w') and B = A 2 /Q. We now assume that the scale over which the correlations of density 

are non-zero is small, in the sense that there is negligible evolution of the density field over the light-crossing time of the correlation scale. 
Writing 

<5(x) = / j^y S k (w) exp(-ik ■ x) (A5) 
and using homogeneity of S to define the power spectrum of the density field P, 

{5 k (w)6* k ,{w')) = (2i:) 3 P(k y w)6 D (k-k') (A6) 
where 8 D is a Dirac delta function. Hence 

(«eff( 9 )K* ff ( 9 )) = B dw / dw' 4>(w)(f>(w')H(w, w', j 9 — 6 [) / dw\ I dw2F(w, Wi)F(w' , w 2 ) 
Jo Jo Jo Jo 

d 3 k 



(2n) 



■ exp[-ik x • (/k(wi) 9 - f K (w 2 ) 9 )] exp[ifc z (wi - w 2 )]P(k, Wi) (Al) 



where k = (k^, k z ) and 

, M^M^ (A8) 

Now we make the standard approximation that the correlation function of 5 is non-zero only if w\ and w 2 are almost equal. Specifically, 
we assume that F and Jk don't vary much over this scale. Thus we set fxiwi) = fn(w 2 ) in the first exponential, and also approximate 
F(w' , w 2 ) ~ F(w' , w\). The w 2 integral then simplifies, since it is approximately 

dw 2 exp(— ik z w 2 ) = 2nS D (k z ) (A9) 

and we are left with 

(«eff( 9 )K* ff ( 9 )) = B / dw dw' 4>(w)<j>(w')H(w,w' ,\ 9 — 9 ) / dwiF(w, w t )F(w', wi) 
Jo Jo Jo 

^iexp[-ik_L ■ f K (m)(^ - 7'))]exp[ifc zWl ]P(|k ± |,^,). (A10) 
(2tt) 2 

An angle integration of P gives 
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Figure Al. Expected weak lensing signal for a AC DM model with erg = 0.9, for the surveys discussed in the main text. The solid line shows equal weighting 
of galaxies; the dotted line shows the effect of employing the proposed multi-colour survey weighting scheme. 



' ^expHkx-MW)(0 - e ))]P{\k±\)= / —P{k, Wl )j Q [kf K {.w l )\e - e \] (Ail) 



J ( 2n ) Jo 2n 

where Jo is the zeroth Bessel function of the first order. Hence we can write the correlation function as 

{Seff( 9 )«eff( )) = B J dw dw' <f>{w)4>{w')H(w, w', I 8 — 6 I) / dwiF(w,wi)F(w' ,wi) 

Jo Jo 
UfjU > >' 

—P{k,wi)Jo[kf K {wi)\0 - 6 \\. (A12) 

Writing this in a form as close to equation E as possible, by reversing the order of integration, and noting that in the weak lensing limit, shear 

1-1 . . 

and convergence have the same statistical properties; (7( 6 )7*( 9 )) — (fteff( )k^ b { 9 )}, we find: 

— + — +' f°° — > — >' f°° hdk — > — »' 

( 7 ( $ ) 7 *( 9 )) = Bj dwiX(w u I 9 - 6 |) J —P(k,w l )J [kf K (w l )\9 - 9 \}. (A13) 

where 

poo poo 

X(wi,9)= / dw dw'<p(w)<p(w')F(w,wi)F(w',wi)H(w,w',9). (A14) 

In the equal-weighted case, H = 1, and X(w,<j)) simplifies to the product of two equal integrals, each independent of <f> and equal to 

g{w)f K (w)/a(w). 

The effect of the weighting is shown in Fig. Al. For surveys such as COMBO-17 with accurate photometric redshifts, the lensing signal 
changes by up to ~ 3%. For the Sloan photometric survey, the effect is ~ 10%; the greatest effect is for RCS, which has the least accurate 
photometric redshifts. 
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